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Abstract. Preliminary femtoscopic results on identical pions from high statistics data set of Au+Au 
collisions at ^.snn — 200 GeV taken during the fourth RHIC run are presented. The measured three- 
dimensional correlation function is studied at low relative momenta using the Gaussian parametriza- 
tion and the Levy stable parametrization. The latter is expected to better describe the data. As the 
results show, both parametrizations underestimate the peak of the measured correlation function 
equally. 
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INTRODUCTION 

Motivation for this study is to check whether a recently proposed parametrization of 
the correlation function using Levy stable source distribution |T] brings a significant 
improvement over the standard Gaussian fit |2]. In addition to this, the non-Gaussian 
source distribution function is also used in most of the models, but the standard method 
of fitting experimental correlation function assumes a Gaussian source 

Possible methods of studying the non-Gaussian effects in the experimental corre- 
lation function include Edgeworth expansion (jSJ and Levy stable source distribution 
parametrization [JJ. 



CORRELATION FUNCTION OF IDENTICAL PIONS 



Event and particle selection criteria. We briefly list the values of the event and 
particle selection criteria used in the present analysis. Detailed description can be found 
in 10]. Events are binned by centrality in five bins corresponding to 0-5%, 5-10%, 
10-20%, 20-30% and 30-80% of the total hadronic Au-i-Au cross-section. In addition 
to the track cuts listed in measured specific ionization of pions is required to be 
farther than ±2a from the Bethe-Bloch theoretical value for electrons. This cut removes 
a contamination due to conversion electrons in the low momentum region. Pairs of 
identical pions are binned by average transverse pair momentum kj = ^\pji + pt2| in 
four bins corresponding to kj e (0.15-0.25,0.25-0.35,0.35-0.45,0.45-0.60) GeV/c 
and the results are presented as a function of the average kj in each of these bins. 



Experimentally, two-particle correlations are studied by constructing the correlation 
function as a ratio 

where A{q) is the measured distribution of the momentum difference q = p\ — pi 
for pairs of particles from the same event and B{q) is the corresponding reference 
distribution for pairs of particles from different events belonging to the same event class 
as analyzed event [2] . The two-particle correlation function at low relative momentum 
q of identical pion pairs is studied using Bertsch-Pratt parametrization iH, @, 0] in the 
longitudinal co-moving system (LCMS) frame, where the relative momentum vector is 
decomposed into the out, side and long components q= {q^.q^, qx) . 



Gaussian parametrization 



Standard method of fitting the two-pion correlation function assumes the Gaus- 
sian source distribution. The correlation function is usually parametrized by a three- 
dimensional Gaussian in ^ 121]. Taking into account a repulsive Coulomb interaction 
between charged identical pions, the measured correlation function dU) is fitted using 
Bowler-Sinyukov procedure ^Q^, 



C{q) = {\-?i)+?iK, 



l+exp(^- 



(2) 



Here, the correlation strength A equals the fraction of pairs originating in the same 
spatio-temporal region relevant for Bose-Einstein correlations, is^ the squared 
Coulomb wave-function integrated over the source with radius 5 fm [4]. are the 
Gaussian source radius parameters defined as the widths of the source emission func- 
tion. Let us note, that only the pairs obeying Bose-Einstein statistics are considered to 
Coulomb interact. For an azimuthally integrated analysis the cross-terms vanish L2J, 



Non-Gaussian parametrization 

Detailed analysis of the shape of the correlation function is important because it 
carries information about the space-time structure of the particle emitting source JUS]- 
Deviations from Gaussian shape can be studied using Edgeworth expansion or Levy 
stable source distribution. 

Edgeworth expansion [3] provides model-independent approach for an analysis of the 
shape of the correlation function. In our previous pion interferometry analysis [4] it was 
shown that the Edgeworth expansion, based on an experimentally preferred Gaussian 
weight function and a complete orthogonal set of even order Hermite polynomials, up to 
6* order is sufficient to describe the data. However, physical interpretation of the higher 
order (4*, 6*) fit parameters is not clear. 



Levy stable source parametrization 



To study possible deviations of the correlation function from a Gaussian shape, we 
followed the method suggested in iQl]. This formalism is relevant for the femtoscopic 
studies of the expanding systems created in heavy-ion collisions, where the scale of the 
fluctuations may be characterized by long tails and asymptotic power-law like behavior. 
The probability distribution of particle emission points then corresponds to a Levy 
stable distribution 111]. Using the Bowler-Sinyukov procedure to determine the repulsive 
Coulomb interaction between identical pions, the two-particle correlation function is 
then characterized by a stretched exponential parametrization. 



C{q) = {\-X) + XK^ 



1 +exp 



a/2 



(3) 



The additional parameter a, when compared to Q, is the Levy index of stability, 
< a < 2. For a = 2 the Gaussian form © is recovered, while for a <2 the correlation 
function becomes more peaked than a Gaussian and develops longer tails. For the 
azimuthally integrated analysis, =0, i 7^ j. 



DISCUSSION OF RESULTS 

Here we present results on two-particle correlations of charged identical pions in Au-i-Au 
collisions at ^^nn = 200 GeV measured in the STAR detector during the fourth RHIC 
run (2003-2004). Data set of 11 x 10^ minimum-bias triggered events is used for this 
analysis. 

The interferometric parameters, correlation strength A and radii Rq, R^ and Ri, are 
obtained by fitting the measured correlation function ^ with the Gaussian parametriza- 
tion Q. Interferometric radius parameters measure the sizes of the homogeneity regions, 
regions from where the particles are emitted with the same average pair momentum kr, 
and their kj dependence contains dynamical information of the pion emitting source jj]. 

Figure [U shows STAR preliminary results on interferometric parameters as functions 
of kj for five centrality bins, where the three-dimensional experimental correlation 
function is fitted using the Gaussian parametrization of the correlation function @. 
Results of the analysis of higher statistics data set of Au-i-Au collisions at ^/¥mn = 
200 GeV are compared with the previously analyzed STAR data 0] for the same system 
at the same energy taken during the second RHIC run (2001-2002). It can be seen that 
the extracted interferometric radii Rq, R^ and Ri are consistent within errors with the 
previous analysis. A small systematic shift in radii can be attributed to the momentum 
resolution correction which is not included in this analysis. Significant difference in A 
is observed in the lowest kj bin. This is explained by an improved purity of the pion 
sample, where the additional cut on particle specific ionization removes contamination 
from the conversion electrons in the low momentum region. 

The Gaussian parametrized fit Q and the Levy stable source parametrized fit Q to 
the measured correlation function, each subtracted from the measured three-dimensional 
correlation function, are projected in the out, side and long coordinates and compared 



in Figure 121 Correlation function is shown for the kj bin 0.25-0.35 GeV/c, for the 
most central collisions and the projections are constrained by the unprojected variables 
qo,qg,qi< 30 MeV/c. It can be seen that contrary to [ 1], Levy stable source distribution 
parametrization Q does not fit the three-dimensional experimental correlation function 
significantly better than the standard Gaussian parametrization ©. Both parametriza- 
tions equally underestimate the peak value of the measured correlation function for 
relative momenta ^o,<?s,^i < 20 MeV/c and underestimate the tail of the correlation 
function. Both effects are mostly visible in the long projection in Figure |2l 

In Figure |3] the interferometric parameters obtained from the Gaussian parametriza- 
tion O are compared to the parameters obtained using Levy parametrization Q of the 
measured three-dimensional correlation function. Interferometric parameters are shown 
as functions of kj for the most central and peripheral collisions. The Levy fit returns sig- 
nificantly larger values of the fit parameters Ro,Rs, R\ and A . However, these parameters 
can not be directly compared to the Gaussian ones representing the interferometric radii, 
because the non-Gaussian parameters do not satisfy the definition of being the widths 
of the source emission function. The large values of the non-Gaussian fit parameters are 
strongly anti-correlated with rather low value of the Levy index of stability a, which 
stays between 1.2 < a < 1.5. 

It can be seen that the Levy stable source distribution parametrization © does not 
bring an advantage in describing the detail shape of the measured correlation function, 
nor in the number of the fitting parameters. Therefore use of the standard Gaussian 
parametrization Q is sufficient and preferred. 

SUMMARY 

The preliminary results on identical pion interferometry using high statistics sample of 
Au-i-Au collisions at ^/^nn = 200 GeV from the STAR experiment at RHIC have been 
presented. 

The results of the Gaussian fit to the measured three-dimensional correlation function 
are consistent within errors with the previously analyzed STAR data. 

It has been shown that in the low relative momentum region the Levy stable source 
parametrization does not fit the experimental correlation function significantly better 
when compared to standard Gaussian parametrization. 

Edgeworth expansion provides the detailed fit to the measured correlation func- 
tion ||4|], but the interpretation of the higher order fit parameters is not clear. 

It seems that to represent the experimental correlation function in Au-i-Au collisions 
at RHIC, the Gaussian parametrization is sufficient. 
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FIGURE 1. Interferometric parameters as functions of and centrality. Results of the Gaussian 
parametrized fit (|2j to the experimental correlation function, Au+Au collisions at y^SNN = 200 GeV. 




FIGURE 2. Gaussian parametrized fit (|2} to data compared to the Levy stable source parametrized 
fit (13 to data, each subtracted from the experimental correlation function. 
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FIGURE 3. Results on interferometric parameters as functions of kj and centrality obtained using 
Gaussian parametrized fit (|3 compared to the parameters of Levy stable source parametrized fit (|3} of the 
measured correlation function. 
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